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Bright solitary waves of atomic Bose-Einstein condensates under rotation 
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We analyse the rotation of bright solitary waves formed of atomic Bose-Einstein condensates with 
attractive atomic interactions. By employing a variational technique and assuming an irrotational 
quadrupolar flow field, we map out the variational solutions in the rotating frame. In particular, 
we show that rotation has a considerable stabilising effect on the system, significantly raising the 
critical threshold for collapse of the bright solitary waves. 
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In recent years, bright solitary waves have been cre- 
ated using ultra-cold atomic Bose-Einstein condensates 

(BECs) [ii,ai. Under attractive atomic interactions, 

these matter waves are self-trapped in the longitudinal 
direction [1, H, @, 0, Q and are closely analogous to the 
classic one-dimensional soliton Q. However, these states 
must be confined in the remaining directions by a waveg- 
uide and can retain three-dimensional effects. The most 
lucid example is the collapse instability: in 3D a homo- 
geneous, untrapped BEG with attractive interactions is 
unstable to collapse JSM while the ID limit is stable to 
collapse. The presence of external trapping stabilises 
the BEG up to a critical atom number (or interaction 
strength) before c ollap se is triggered, as demonstrated 
experimentally [O, [l^- The collapse instability has lim- 
ited bright solitary wave (BSWs) experiments to only a 
few thousand atoms per BSW Q,[3,ti. The BSW solu- 
tions, and their critical points, have been studied theo- 
retically 11, i, 0, l, 

with variational approaches shown 
to give very good predictions. The collisions of BSWs, 
which show intriguing behaviour and may have applica- 
tions in interferometry, are also prone to collapse instabil- 
ities 0, H, [13, [IE [Hi • As such, it is pertinent to consider 
approaches to suppress collapse in attractive BEGs. Pe- 
riodic modulation of the interaction strength, made pos- 
sible by employing a Feshbach resonance, is predicted 
to partially stabilise against collapse [l^ , although when 
the average interaction is attractive, collapse is inevitable 
[l7| . The presence of a quantized vortex is also predicted 
to raise the threshold for collapse [l^, [l^ , although the 
presence of a vortex in an attractive condensate under 
harmonic trapping is not energetically stable [H, US] • 

Due to the superfluid nature of BEGs, it is intriguing 
to study their response to rotation. The rotation of repul- 
sive BEGs has been considered extensively both experi- 
mentally and theoretically (see [2l| for a review). One 
method is to mechanically rotate the BEG in an elliptical 
trap [55, '23', '23| formed by time-dependent laser or mag- 
netic fields. At low rotation frequency Q, the condensate 
remains irrotational and vortex-free. According to a hy- 
drodynamical model, the BEG can access a family of ro- 
tating stationary solutions, characterised by a quadrupo- 
lar irrotational flow pattern [H IH, 113, [ll, [ll, ^ and 
confirmed experimentally [1^. At a critical rotation fre- 



quency, which coincides with when these irrotational so- 
lutions become unstable [1^, [s^] , vortices are nucleated 
and form a vortex lattice. In the context of attractive 
BEGs, theoretical work has shown that a centre-of-mass 
mode is favoured under rotation rather than the oc- 
curence of vortices [10, [3l|. Under harmonic trapping, 
this mode becomes excited when Q exceeds the trap fre- 
quency. 

In this work we consider the rotation of bright solitary 
matter waves in an elliptical trap about the longitudinal 
axis. We employ a variational technique based on assum- 
ing an ansatz for the BSW profile which incorporates a 
quadrupolar irrotational flow pattern. By deriving the 
variational energy of the system, we obtain the BSW so- 
lutions and analyse their response to rotation. In partic- 
ular, we find that rotation of the BSW can significantly 
increase the critical point for collapse. 

We consider a BEG confined by an atomic "waveguide" 
potential, under rotation at frequency f2 about the longi- 
tudinal axis. In the limit of zero temperature, the BEG 
can be described by a mean-field "wavefunction" ^ (r, t ) 
which satisfies the Gross-Pitaevskii equation (GPE) [2]| . 
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where m is the atomic mass and g = ATrh^as/m parame- 
terizes the atomic interactions, with as being the s-wave 
scattering length. The term accounts for frame ro- 
tation, where — —ih{x-^ — y-§^) is the z-component 
angular momentum operator. We assume that the con- 
fining potential ^(r) is harmonic with the form. 
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where oj^ is the average trap frequency in the x-y plane, 
e is the trap ellipticity in the x-y plane, and A = uj^/^r 
is the trap ratio that determines the axial trap strength. 

In order to study the BSW solutions we employ a vari- 
ational technique. This involves assuming a BSW ansatz 
and minimising its energy to obtain the variational solu- 
tions. This technique has been employed for non-rotating 
BSWs and trapped attractive BEGs, and has been shown 
to have give very good agreement with the full solution 
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of the GPE. When the axial trapping is weak (A ^ 1) 
we will employ the sech ansatz for the BSW, given by, 
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where N is the atom number. Ignoring the e*"^'' term, 
this ansatz is identical to that used in Refs. [1, 0, d] , with 
the sech axial profile appropriate because it is the form of 
the ID soliton solution [9^. The term e'"^^ introduces a 
quadrupolar flow pattern to the BEG, with a determining 
the amplitude. This flow pattern preserves irrotationality 
and has been very successful in modelling vortcx-frec re- 
pulsive condensates under rotation [23, 26, 27, 28, 29, 30]. 
The total energy of the system is defined by. 
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Insertion of the ansatz into Eq. ^ gives the energy iJs, 
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We can reduce the number of variables in Eq. ^ as 
follows. Under the hydrodynamical interpretation, the 
mean- field wavefunction can be expressed as ^' (r, t) = 
^yn{r,t) exp[i(^(r, t)] where n(r, t) and 4>{r, t) are the con- 
densate density and phase, respectively. Furthermore, 
v(r, i) = {h / m)W (f>{r , t) is the "fluid" velocity. Insert- 
ing this into the GPE and equating imaginary parts, one 
derives a continuity equation given by. 



dn 
'dt 



+ V • (n [v- rj X r]) = 0. 
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Our irrotational phase distribution 0(r) = axy corre- 
sponds to a velocity distribution v(r) = fia{y\ + xj)/m. 
Inserting this into the continuity equation and setting 
dn/dt — 0, we find that stationary solutions satisfy. 
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Note that a can be positive or negative, resulting in two 
"branches" of solutions. For a > the BSW is wider 
in the a;-direction than in the ^/-direction, and vice versa 
for a < 0. We can thus eliminate a from Eq. ([5]). For 
simplicity we employ rescaled variables (in terms of the 
transverse harmonic oscillator) = 'x/or, ly = ly/ar, 
'j^ = Iz/ttr, f2 = ri/wr and es — Es/{NhuJr), where 



flj. = y^h/ mujr is the radial harmonic oscillator length. 
Furthermore, we introduce the dimensionless interaction 
parameter k — N\as\/ar. The ansatz energy becomes, 
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Recall that Eq. ^ is vahd for A < 1. Under tight 
axial trapping A 3> 0, this direction is dominated by the 
trap rather than the interactions and it is more appro- 
priate to consider a gaussian ansatz, 
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The rescaled energy for the gaussian ansatz is then. 
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Note that, for the regimes of interest, the sech and gaus- 
sian ansatzgive similar results, typically differing by less 
than 10% We can now map out the 3D energy land- 
scapes of the rotating BSWs as a function of the length- 
scales 7a;, 7y and 7z. Variational BSW solutions exist 
where there is a local energy minimum in the landscape, 
and are obtained by a numerical search algorithm. The 
local energy minimum has widths 7^, 7]^ and 7^, energy 
So, and quadrupolar flow amplitude ap. 

We first revisit the A = BSW solutions in the absence 
of rotation and ellipticity, as studied previously using the 



(a) 




0.5 1.0 

Yr 



1.5 




FIG. 1: Energy landscape of the non-rotating and non- 
elliptical system for A = according to Eq. (|8]). (a) Stable 
regime k — 0.4 < kc featuring a local energy minimum, i.e. 
the BSW solution, (b) Unstable regime k — 0.8 > kc, where 
the whole parameter space is unstable to collapse. White 
contours highlight the shape of the landscapes. Since Eq. (|8]) 
is cylindrically symmetric in this case, we introduce a radial 
lengthscale jr = \/ (7I -I- 7h)/2. 
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FIG. 2: Properties of e = BSW solutions as a function of fl 
for A = (solid line) and the axially-trapped cases of A = 1 
(dashed lines) and A = 5 (dotted lines) . (a) Quadrupolar flow 
amplitude ao for interaction parameter k = 0.4. (b) BSW 
energy eo for k — 0.4. (c) Axial lengthscale 7^ for k — 0.4. 
(d) Critical interaction parameter for collapse kc. Note that 
for A = (A > 0) we employ the sech (gaussian) ansatz. 



a = limits of Eq. dH]) and Eq. ^ ^. Here the 

energy landscape is cylindrically symmetric (e = 0) and 

so we introduce a radial lengthscale 7^ = ^^(7^ + 7^)/2. 

A typical energy landscape, according to Eq. ([5]), for 
a stable BSW solution is presented in Fig. 1(a), corre- 
sponding to fc = 0.4. At the origin the interaction term in 
Eq. ([8]) diverges to negative values and is a region of col- 
lapse of the BSW. However, there exists a local energy 
minimum which represents the self-trapped BSW solu- 
tion. A typical unstable energy landscape is shown in 
Fig. 1(b) for k — 0.8. No local energy minimum exists, 
and the whole parameter space is unstable to collapse. 
From this the critical interaction strength for collapse is 
determined to be fcc = 0.76, in good agreement with full 
solution of the GPE which gives « 0.68 

We will now consider the effect of rotation. For sim- 
plicity we first assume e = 0. Figure 2(a)-(c) shows how 
the key parameters vary as rotation is introduced for a 
fixed interaction parameter k — 0.4. For > the sym- 
metry between 73, and "/y is broken and the solutions 
have non-zero ao (Fig. 2(a)). We see the formation of 
two branches of Uq. Due to the trap symmetry in the 
x-y plane, the branches are symmetric about the ao = 
axis, with the upper branch being elongated in the x- 
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FIG. 3: Properties of the rotating A = BSWs for ellip- 
tical traps of e = 0.1 (dashed lines) and 0.2 (dotted lines). 
Black (grey) lines indicate upper (lower) branch solutions, (a) 
Quadrupolar flow amplitude qq for k = 0.4. (b) BSW energies 
£0, rescaled by their non-rotating values of eo{^ = 0) = 0.972 
for A = 0, 1.315 for A = 1 and 3.01 for A = 5. (c) Axial 
lengthscales rescaled by 7^(0 = 0) = 2.36 for A = 0, 0.902 
for A = 1 and 0.416 for A = 5. (d) Critical interaction pa- 
rameter for collapse fcc of the upper branch (lowest energy) 
solutions. 



direction and the lower branch being elongated in the 
7/-direction. As ft increases, so too does the magnitude 
of ao, implying a spreading of the BSW in the x-y plane. 
This is because of the growth of an outward centrifugal 
force. As fl approaches ujr, ao diverges to ±00. This 
is because, at n — ujr, the centrifugal force exactly bal- 
ances the trapping potential, and the BEG is untrapped 
in the x-y plane. Since the BEG centre-of-mass becomes 
dynamically unstable, this is termed the centre-of-mass 
instability [33] • The BSW energy (Fig. 2(b)) decreases 
towards zero as fl —^ u>r as a result of the reduced den- 
sity. The axial lengthscale (Fig. 2(c)) grows with since 
the radial spreading dilutes the interaction strength and 
forms a less tightly bound BSW. 

We have isolated the critical interaction parameter for 
collapse kc as a function of Q with the results shown in 
Fig. 2(d). The most striking feature is that kc dramat- 
ically increases as flz approaches ujr- This is directly 
associated with the radial spreading and reduced density 
of the rotating solutions. Specifically, for U/ujr = 0.9, kc 
is approximately 50% larger than its non-rotating value 
while for Q/ujr = 0.97, kc is approximately twice as large. 

In Fig. 2 we also consider the presence of significant 
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axial trapping A = 1 (dashed line) and 5 (dotted line), 
for which we employ the gaussian ansatz. We see simi- 
lar qualitative behaviour to the A = case: a divergent 
growth of ao (Fig. 2(a)) and decrease in Eq (Fig. 2(b)). 
However, the magnitudes are consistently less than the 
corresponding A = results. The axial lengthscales (Fig. 
2(c)) show little variation with since this is now dom- 
inated by the external axial trapping. The critical point 
for collapse kc grows with 17, but at a slower rate than for 
A = 0. Note that the presence of axial trapping reduces 
fcc, as observed elsewhere [sl. [isj. 

The e = limit is somewhat unphysical since no torque 
is actually apphed to the BEC. We now consider the more 
realistic case of finite trap ellipticity. The results for A — 
BSWs under e = 0.1 and 0.2 are presented in Fig. 
3. The finite ellipticity breaks the symmetry in the x-y 
plane and therefore in the branches of ao. The upper 
branch solutions are elongated in the i-direction, which 
has trap frequency uj^ = ^(1 — e)LOr. The centrc-of-mass 
instability then occurs when 51 = ■\/(l — €)ujr, which is 
why the divergence in ao shifts to lower il. Conversely, 
the lower branch solutions diverge at SI = ^/l + eujr and 
so become shifted towards larger SI. 

Consideration of the energy (Fig. 3(b)) shows that the 
upper branch solutions have lower energy. That is, it is 



lower energy for the condensate to be elongated in x than 
y, since the trap is weaker in this direction. Consequently, 
we expect that only the upper branch solutions would 
ever be observed. Apart from the shift in the asymptotes 
introduced by the finite ellipticity, the behaviour of ao, eo 
and A^ is qualitatively and quantitatively similar to the 
e = case. In Fig. 2(d) we plot the critical interaction 
parameter for collapse kc of the upper branch solutions. 
Again, we see similar behaviour to the e = case, with a 
dramatic in increase in k^ as Q, approaches ^Jl — eto^- 

In this work we have employed a variational technique 
to study bright solitary matter-wave solutions under ro- 
tation in elliptical traps. This is made possible by incor- 
porating an irrotational quadrupolar flow pattern into 
the variational ansatz. Importantly, the BSW becomes 
more stable to the collapse instability when under rota- 
tion. This is most pronounced when the rotation fre- 
quency SI is slightly less than the minimum effective trap 
frequency — eujr, e.g., for fl — 0.9^/1 — etUr, the crit- 
ical interaction parameter for collapse kc increases by 
approximately 50% over the non-rotating value. This 
implies that the BSW can support the same increase in 
number of atoms, before collapse. 
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